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ABSTRACT 

We investigate here the effects of plasma instabilities driven by rapid pair cascades, which arise in the 
environment of GRB sources as a result of back-scattering of a seed fraction of the original spectrum. The 
injection of e pairs induces strong streaming motions in the ambient medium. One therefore expects the 
pair-enriched medium ahead of the forward shock to be strongly sheared on length scales comparable to the ra- 
diation front thickness. Using three-dimensional particle-in-cell simulations, we show that plasma instabilities 
driven by these streaming e pairs are responsible for the excitation of near-equipartition, turbulent magnetic 
fields. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective 
coupling between e pairs and ions, and may help explain the origin of large upstream fields in GRB shocks. 
Subject headings: gamma rays: bursts - instabilities - magnetic fields - plasmas - shock waves 



1. INTRODUCTION 

More than three decades ago, it was pointed out that 7-rays produced in suffic iently lumi nous and compact astrophysical 
sources would create e pairs by collisions with lower energy photons: 77 — > e + e~ (Jell ey|l966[ ). This mechanism both depletes 
the escaping radiation and also changes the composition and properties of the radiating gas through the injection of new particles. 
An approximate condition for such pair creation to become significant is that a sizable fraction of the radiation from the object 
be emitted above the electron mass energy, m e c 2 = 511 keV, and that the compactness of the source, 

Lctt 

l=j (1) 

exceeds 1, where L is the total luminosity, and r/ is the characteristic source dimension. When I ~ 1, a photon of energy 
e = hv/m e c 2 ~ 1 has an optical depth of unity for creating an e pair. 

As an illustration consider a spherical source with a spectrum emitting a power L in each decade of frequencies, I would be 
> 1 for 7-rays of energy e if the radius of the source satisfies 

n < 10 17 e( — ^ _, cm. (2) 



v 10 49 ergs 

The relevant values of L range from 10 50 - 10 52 erg s" 1 for GRBs. In addition, the short spikes observed in the high energy 
light curves suggest that GRBs dissipate a significant fraction of their energy at r <C r;, so that the source is indeed so small that 
equation (|2| is easily satisfied. This argument is, however, only applicable to 7-ray photons emitted isotropically and is alleviated 
when the radiating source itself expand at a relativistic speed ( |Piran|1999] l. In this case, the photons are beamed into a narrow 
angle 6 ~ 1/T along the direction of motion, and, as a result, the threshold energy for pair production within the beamed is 
increased to e ~ V. 

The non-thermal spectrum of GRB sources is therefore thought to arise in shocks which develop beyond the radius at which 
the relativistic fireball has become optically thin to 77 collisions ( |Piran 1 1 999) . However, the observed spectra are hard, with a 
significant fraction of the energy above the 77 — > e + e~ formation energy threshold, and a high compactness parameter can result 



in new pairs being formed outside the originally optically thin shocks responsible for the primary radiation ( |Thompson & Madau 
2000 ). Radiation scattered by the external medium, as the collimated 7-ray front propagates through the ambient medium, would 
be decollimated, and, as long as / > 1, absorbed by the primary beam. An e pair cascade can then be produced as photons 
are back- scattered by the newly formed e^ pairs and interact with other incoming see d photons (|Thompson & Madau ||2000 



|Beloborodov|2002| |2005] Meszaro s^d.|2001t|Ramirez-Ruiz et al. 2002 ; |Li et al.|2003||Kumar & Panaitescu|2004) . 

In this paper, we consider the plasma instabilities generated by rapid pair creation in GRBs. The injection of e^- pairs 
induces strong streaming motions in the ambient medium ahead of the forward shock. This sheared flow will be Weibel-like 



unstable (Medvedev & Loeb 1999), and if there is time before the shock hits, the resulting plasma instabilities will generate sub 
equipartition quasi-static long-lived magnetic fields on the collisionless temporal and spatial scales across the e ± pair-enriched 
region. This is studied in ^4] using three-dimensional kinetic simulations of monoenergetic and broadband pair plasma shells 
interpenetrating an unmagnetized medium. The importance of the electromagnetic filamentation instability in providing an 
effective coupling between e ± pairs and ions is investigated in ^2] while a brief description of the numerical methods and the 
initial models is giving in ^3] The implications for the origin of the upstream magnetic field, in particular in the context of 
constraints imposed by observations of GRB afterglows, are discussed in ^5] 
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2. PAIR LOADING AND MAGNETIC FIELD GENERATION 

Given a certain external baryon density n p at a radius r outside the shocks producing the GRB primordial spectrum, the initial 
Thomson scattering optical depth is r ~ n p ojr and a fraction r of the primordial photons will be scattered back, initiating a pair 
e cascade. Since the photon flux drops as r" 2 , for a uniform (or decreasing) external ion density most of the scattering occurs 
between r and r/2, and the scattering and pair formation may be approximated as a local phenomenon. 

Consider an initial input GRB radiation spectrum of the form F(e) = Fp(e/ep)~ a for e > e«, where ~ 0.2- 1.0 is the break 
energy above which the spectral index a ~ 1 — 2 (for the present purposes the exact low energy slope is unimportant). Radiation 
scattered by the external medium is therefore decollimated, and then absorbed by the primary beam.The impact of such process 
strongly depends on how many photons each electron is able to scatter ( |Beloborodov|2002] l 

A r ~4x 10 8 



(l0^) ( 1051 erg s-' ) Cm<A ' (3) 



10 51 erg s 

where \j = Xjin^aj) w c/iFpaj) is the electron mean free path and n 7 is the photon density. The photons forced out from 
the beam by the electrons can then produce pairs as they interact with incoming photons (Fig. Ill, so that a large number of 
scatterings implies a large number of pairs created per ambient electron. 
At some distance G3 < A from the leading edge of the radiation fr ont, the numbe r of ph otons scattered by one ambient electron 



is ~ B3/A7- and a fraction ~ B3/A 77 of these photons are absorbed peloborodov (2002J). One e pair per ambient electron is 



1987) and e th is the 



injected when G5 = (A r A 77 )'/ 2 « \ T e" h /2 /[Il(a)] 1/2 > X T , where n(a) = 2" Q (7/12)(1 +aT 5 ^ ( |Svensson| 
photon threshold energy for e ± formation. For 1.5 < a < 2, one has 15 A7- < 05/ < 25Aj-, so that pair creation substantially lags 
behind electron scattering ( Beloborod ov|2002 1 



Most of the momentum deposited through this process involves the side-scattering of very soft photons, which collide with 
hard 7-rays to produce energetic (and almost radially moving) pairs. A photon of energy t,< 1 that is side-scattered through 
an angle ~ 9 r creates a pair if it collides with another photon with energy exceeding e t h ~ 4(6* 2 e r )~ I . The injected pair will be 
relativistic with j± ~ e t h > 1. The distribution of injected pairs is therefore directly determined by the high energy spectral index 
a. This motivates our study in ^4]of radially streaming, relativistic pair plasma shells with a broadband kinetic energy distribution 
interpenetrating an unmagnetized medium. 

This pair-dominated plasma, as long as its density n± < n p (m p /2m e ), is initially h eld back by the inertia of its constituent 
ions, provided that the pairs remain coupled to the baryons (Thompson & Madau 2000; Beloborodov 2002). The latter is 



likely to be the case in the presence of weak magnetic fields (Thomps on & Madau ||2000 i. In the absence of coupling, the 
pair density would not exponentiate, mainly due to the (1-/3) term in the scattering cross section, where (3 = v/c. Instabilities 
caused by pair streaming relative to the medium at rest, as we argue in ^4] are able to generate long-lived magnetic fields on the 
collisionless temporal and spatial scales, which are modest multiples of the electron plasma frequency, uj e = (Aire 2 n e /m e ) x l 1 , and 
the collisionless skin depth, 

c . _5 / n 



\ e = — ~5xl0 3 (— ^-1 ' /2 <A r <A 77 . (4) 
uj e VI cm^ / 

Here n e is the electron number density. Plasma instabilities therefore evolve faster than the time between successive scatterings 
and much before the scattered photons are absorbed by the primary radiation (Beloborodov 2002). In the presence of transverse 



magnetic field B the pairs gyrate around field lines on the Larmor time, ujg 1 = m e c/(Be). The net momentum of the e ± pairs is 
thus efficiently communicated to the medium (Frederiksen et al. 2004; Hededal & Nishikawa 2005). Magnetic coupling may 
dominate if u>b > w c , which requires B 2 /4ir > n e m e c L ( |Beloborodov|2002| . In the sections that follows we present a quantitative 
discussion of the effects of plasma instabilities driven by rapid e pair injection. 

3. SIMULATION MODEL 

Here we illustrate the main features of the collision of an e ± pair plasma shell into an unmagnetized medium, initially at 
rest, using a modified version of the PIC code TRISTAN - first developed by Buneman ( 1993) and most recently updated by 
Nishikawa et al. ( |2005| [2006). The simulations were performed on a 85 x 85 x 640 grid (the axes are labeled as x, y, z) with a 



total of 3.8 x 10 particles with periodic ([x,y] plane) and radiative (z direction) boundary conditions. In physical units, the box 
size is 8.9 x 8.9 x 66.7 (c/aj e ) 3 , and the simulations ran for 6Qoj~ 1 . 

In the simulations, a quasi-neutral plasma shell, consisting of e ± pairs and moving with a bulk momentum u- = 7oV z /c along 
the z direction, penetrates an ambient plasma initially at rest. Here 70 is the initial Lorentz factor of the pairs and v z is the bulk 
velocity of the shell along z. e pairs are continuously injected at z = 2.6A e = 25A, where A = X e /9.6 is the grid size. The 
ion-to-electron mass ratio of the ambient plasma is set to m, /m e = 20, with both plasma populations having a thermal spread with 
an electron rms velocity v^/c = 0.1. For completeness, a purely pair ambient plasma (i.e., a mass ratio mjni e of 1) is also 
studied, with a thermal electron velocity \%/c = 0.1. The shell and the ambient plasma have a density ratio of 0.75. Two different 
bulk Lorentz factor configurations for the injected (cold) e ± pairs are considered: a monoenergetic (u z = 12.5, v±/c = 0.01) 
and a broadband distribution (3.0 < u z < 30.0, v t h/c = 0.01). Both distributions have similar kinetic energy contents and plasma 
temperatures (Fig.[2j). 

4. RESULTS AND INTERPRETATION 

Much of our effort in this section is dedicated to investigate the ability of e pair loading to alter the physical parameters 
characterizing the ambient medium. Some of the questions at the forefront of attention include the level of the (electro) magnetic 
field generated via plasma instabilities, as well as the saturated state of the particles and fields. We address all of these issues 
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here. We consider both e pair and electron-ion plasmas as possible ambient media. In the dimensional units used, the results 
for e pair ambient plasmas are equally applicable to interpenetrating proton/antiproton plasma shells. 

4.1. e ± Pair Ambient Plasma 

As pointed out in |Silva et aL| ( |2003) >, encountering the medium at rest, the incoming pairs are rapidly deflected by field 
fluctuations growing because of the two-stream instability (Me dvedev & Loeb|1999 Pruet et al.|2001| Gruzinov|2001[ ). A large 



number of oppositely directed current filaments are generated, as illustrated in Fig. |3| which in turn will generate inhomogeneities 
in the magnetic field (predominantly in the x-y plane). This configuration is unstable because opposite currents repel each other, 
whereas like currents are attracted to each other and tend to coalesce and form larger current filaments. During the linear stage of 
the instability, the rapid generation of a large number of randomly distributed current and magnetic filaments is observed (with 
a correlation length Ay w 2 x l A ^fxh(c/uj e ) = 1.2X e = 10A) , which in turn is accompanied by a expeditious production of a strong 
magnetic field (Fig.[3]l. The magnetic field energy density reaches 0.9% of the initial total kinetic energy (e p ) in the broadband 
case and 1.2% in the monoenergetic one (i.e., e# = J {B 2 dV /Sir)/ e p =0. 009 and 0.012, respectively). 

As the instability enters the saturation phase, these initially randomly oriented filaments begin to interact with each other and 
merge in a race where larger electron channels consume smaller neighboring channels. In this manner, the transverse magnetic 
field grows in strength. The instability then saturates and the energy in the magnetic field decays. The strong decrease in the 
magnetic field energy density is associated with a topological adjustment in the structure of currents and fields (Silva etaL] 



2003 ), and is mostly linked to a reduction of the field's volume filling factor. These distinctive features are present in both the 
monoenergetic and broadband cases (Fig. [3]). However, there are some clear contrasts. For a plasma shell with a broadband 
distribution of initial Lorentz factors, transverse energy spreading happens over a variety of timescales. This can be understood 
as fol lows. The instability aris es from the free streaming of particles, with a corresponding linear growth rate scaling with 



-1/2 

7 (Medvedev & Loeb 1999). Our experiments show that this does indeed happen; the continuous injection of a broadband 



distribution of momentum leads to a diversity of linear growth rates. As the field amplitude grows, the transverse deflection 
of particles increases (with decreasing 70), and thus free streaming across the field is suppressed preferentially for the mildly 
relativistic pairs. The magnetic field energy grows in the early stages by slowing down the pair plasma shell. 

Satu ration is then achi eved by the combination of transverse energy spreading and the generation of near-equipartition magnetic 
fields (Silva et al.|2003) >. In the broadband case, transverse energy s preading happens more swiftly for mildly relativistic pairs, 
and, as a consequence, B sat is decreased. This follows directly from (Me dvedev & Loeb|l999| l the scaling of B sdt with v^ 1 (or in 
this case, v^ 1 )- As can be clearly seen in Fig. |4j the ambient e ± pairs are more effectively heated in the broadband case. After 
saturation, which takes place earlier for mildly relativistic pairs, the energy stored in the magnetic field is transferred back to the 
plasma particles, leading to strong heating and the generation of a high-energy tail in the distribution (Fig. El. The e ± pairs are 
nonetheless expected to thermalize, given sufficient space and time (Spitkovs y |2006| l. 

It is not surprising that the linear growth rate saturates at comparable timescales (Fig. |2j, since in the broadband distribution, 
the energy averaged Lorentz factor, 70, is not that different from 70 in the monoenergetic case. Fig. [3] also shows a transient, 
electromagnetic precursor produced as the first injected e pairs stream through the ambient medium unhampered. The saturation 
level of the magnetic field in the precursor region is expected to occur over longer distances. This is because the linear growth 
rate of the electromagnetic filamentation instability in the precursor region is greatly decreased as a result of e ± pairs being 
continuously depleted from the shell's front (Spit kovsy |2006 1. 



4.2. Electron-Ion Ambient Plasma 

Encountering the medium at rest, the incoming e ± pairs are rapidly deflected by field fluctuations growing because of the 
two-stream instability. The initial perturbations grow nonlinear as the deflected e pairs collect into first caustic surfaces and 
then current channels. The resultant cylindrical magnetic fields cause mutual attraction between currents forcing like currents 
to approach each other and merge. As a result, the magnetic field grows in strength. This continues until the fields grow strong 
enough to deflect the much heavier ions (Frederi ksen et al.|2004) . Fig.|5]shows that this happens over ~ 20 electron skin depths 
from around z = 300. As illustrated in Fig. [5] the ions stay clearly separated in phase space and are only slowly heated. In the 
presence of ions, the incoming e ± pairs will drive higher levels of saturated B-field (Fig.pl, by a factor of (mj/mj 1 / 2 , albeit on 
a longer timescale: the magnetic field energy density in the monoenergetic case reaches 6.1% of the initial total kinetic energy in 
the electron-ion ambient plasma case and 1.2% in the e pair one. This is due to the massive ion bulk momentum constituting a 
vast energy reservoir for particle heating. 

In the case of an electron-ion ambient plasma, the ion channels are subjected to a similar growth mechanism as the positrons, but 
at a slower rate. When ion channels grow sufficiently powerful, they begin to experience Debye shielding by the electrons, which 
by then have been significantly heated by scattering on the increasingly larger electromagnetic field structures. The temporal 
development of the pair and ion channels is illustrated in Fig. [6] The large random velocities of the electron population allow 
the concentrated ion channels to keep sustaining strong magnetic fields. 

From the requirement that the total plasma momentum should be conserved, the electromagnetic field produced by the two- 
stream instability acquires part of the longitudinal momentum lost by the counterstreaming populations. Consequently, the 
magnetic field energy grows in the early stage by slowing down the e ± pair plasma, as illustrated in Fig. [2] Our experi ments 
show that this also happens for e ± pair shells injected into an electron-ion ambient plasma (Fig|7j. Although as argued in §4.1' 
there are some clear differences in the evolution of the associated absorption of momentum between the broadband and 
monoenergetic cases. Figs. [8] and |9| elucidate some of these differences. 

As can be clearly seen in Fig. pf transverse energy spreading occurs faster in the broadband case, and, as a consequence, the 
ambient ion beam is heated more promptly. The total magnetic energy then grows as the ion channels merge. The magnetic field 



the 
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associated with these currents has also a filamentary structure, as seen in Fig. [9] The magnetic energy scales with the square of 
the electric current, which in turn grows in inverse proportion to the number of current channels. The net result is that the mean 
magnetic energy increases accordingly: the magnetic field energy density reaches 6.1% (4.7%) of the initial total kinetic energy 
in the monoenergetic (broadband) case. 

The ambient ions retain distinct bulk speeds in shielded ion channels and thermalize much more slowly (Fig. [8]l. This is 
expected to continue as long as a surplus of bulk relative momentum remains in the counterstreaming plasmas. Fig. [STshows the 
extension of the ion currents and the corresponding magnetic field energy density. After the linear stage, the instability saturates 
and the energy in the magnetic field decays, reflecting essentially the subsequent decrease in the field's volume filling factor. 
In the monoenergetic case, the magnetic energy density drastically drops after about ~ 30^"' . The decrease in the magnetic 
field energy is, however, slower in the broadband shell. This is associated with the continuous transverse energy spreading of 
progressively faster moving e ± pairs. The transverse energy spreading is less efficient for the faster moving component and, as 
a result, a large number of oppositely directed current filaments continue to be generated even after the bulk of the flow enters 
the saturation phase. In the monoenergetic case, on the other hand, the magnetic energy density drastically drops after about 
~ 3Quj~ 1 . These fields, however, maintain a strong saturated level for at least the duration of the simulations. The final magnetic 
energy density level is still quite high, up to 2.75% (2.57% ) of the initial total kinetic energy in the monoenergetic (broadband) 
case. 

This is in contrast to the results reported by Frederiksen et al. ( 2004| ) for relativistically counterstreaming electron-ion plasma 



collisions, where the magnetic field energy density continues to grow throughout their experiment, which lasts for about 480WJ 1 . 
In these experiments, both counterstreaming plasmas are, however, composed of electron and ions, which in turn significantly 
increases the relative bulk momentum. It is thus clear that the symmetry associated with the counterstreaming electron-ion 
populations is essential for the longevity of the ion current channels. 

5. CONCLUSIONS 

We present self-consistent three-dimensional simulations of the fields developed by the electromagnetic filamentation instabil- 
ity as a consequence of e ± pair injection, which arise naturally in the environment of GRB sources as a result of back-scattering. 
Our results demonstrate that even in an initially unmagnetized scattering plasma, a small-scale, fluctuating, predominantly 
transversal, and near-equipartition magnetic field is unavoidably generated. These fields maintain a strong saturated level on 
timescales much longer than X e at least for the duration of the simulations ~ 60(c/A e ) ~ 0.1(c/Ar). The e ± pairs are effectively 
scattered with the magnetic field, thus effectively communicated their momentum to the scattering medium initially at rest. Our 
results indicate that the fields necessary to ensure that e pairs remain coupled to the medium can be easily created via plasma 
instabilities. The next required step is to increase the ion-to-electron mass ratio in the scattering medium in order to determine 
the spatial spread and character of the particle coupling. 

A question which has remained largely unanswered so far, is what determines the characteristic strength of the upstream 
magnetic field in afterglow shocks, which is inferred to extend to tenths, or even a tens, of mG ( jLi & Waxman|2006] l. This is 
of course large in comparison with the magnetic fields presumed to be present in the interstellar medium, which are measured 
in fj, G. A sufficiently strong magnetic field may be present in a wind ejected by the GRB progenitor but only if its is highly 
magnetized, which is not thought to be the case for the widely favored Wolf-Rayet stars ( |Li & Waxman|2 006). The injection 
of e pairs induces strong streaming motions in the ambient medium ahead of the forward shock. This sheared flow, as we 
demonstrated here, will be Weibel-like unstable and, if e ± loading continues, well after the shock hits the e ± pair enriched 
mediu m, the resulting upstream flow will be strongly magnetized. Since one expects the radiation front to lead the forward 
shock ( |Beloborodov|2002 | |2005| by a small distance, R/(4T 2 ), one expects the instability to be able to generate subequipartition 
magnetic fields, B ~ 0.\(n p /\ cm -3 ) 1 / 2 G, in the scattering medium just before it is shocked. In t his case, the constr aint on the 
upstream field strength of B 0.05(n p /T cm" 3 ) 5 / 8 mG, imposed by X-ray afterglow observations (Li & Waxman 2006) is easily 
satisfied. 
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FIG. 1 . — Schematic plot of pair cascades triggered by the back-scattering of seed 7-ray photons on the external medium. 
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FIG. 2. — e^ pair momentum distribution functions at t = 20.8w~' (solid curves) and t = 59.8u~' (dotted curves). The blue (red) curves are for a monoenergetic 
(broadband) initial momentum distribution function. The top (bottom) panel show the distribution functions for particles moving parallel (perpendicular) to the 
injected e^ pair plasma. 
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FIG. 3. — Growth of the two-stream instability at time t = 59. 8^' for the two different momentum distribution function shown in Fig.ol Here we show the 
average transverse magnetic field amplitude (in simulation units) in the x-y plane as a function of z (solid curves). The various panels show the (x,y) components 
of the magnetic field and the z component of the current density (J z ) in different subsections of the computational box. The top (bottom) panels are for a 
simulation in which the initial momentum distribution function magnetic field is monoenergetic (broadband). Color bar gives the amplitude of J z = [7 ma x > — /max] 
in simulation units for : /max = 20.04 (25.57) for z=150, /max = 99.5 (67.12) for z=350, and / ma x = 35.41 (29.93) for z=520. The expected correlation length, A/, 
of the randomly distributed current and magnetic filaments is plotted in the z=350 subsections. 
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FIG. 4. — Current filamentation and magnetic field energy density after saturation: t = 59. 8u>~ . The amplitude of the average (in the x-y plane) transverse 
magnetic field, B± is plotted (in simulation units) as a function of z together with the average z component of the current density, 7-, for both injected (top) and 
ambient medium (bottom) pairs . The a (b) panels are for a simulation in which the initial momentum distribution function is monoenergetic (broadband). 
Color bar gives the amplitude of J z = [7 m ax, — /max] an d — [0,Z? m ax] in simulation units: Jmax — 18.5 (17.3) for the injected e ± pairs, 7 max = 15.9 (13.9) for the 
ambient pairs and B max = 1.8 (1.2). 
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FIG. 5. — Longitudinal heating and acceleration, illustrated by changes in «n for both injected (monoenergetic) pairs (red) and ambient plasma (blue). The 

bottom (top) panels are for a simulation in which the ambient medium is composed of electrons and ions (e pairs). In the top panel, the ambient positrons are 
initially at rest but are strongly accelerated by the jet. In the bottom panel, the ions, being heavier than the pairs, remain clearly separated in phase space and are 
only slowly heated. Also shown are the average transverse magnetic field amplitude (in arbitrary units) in the x-y plane as a function of z (solid curves). 
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FIG. 6. — Longitudinal current densities in the counterstreaming plasmas. The various panels show the z component of the current density, J z , in different 
subsections of the computational box for pairs injected with a monoenergetic distribution. The small insets show the ion current in the same plane. Color bar 
gives the amplitude of J z = [7 maX! — J max ] in simulation units for pairs (ions) : J m!lx = 20.04 (0.1) for z=150 and 7 max = 173.3 (8.1) for z=350. The expected 
correlation length, Xj, of the randomly distributed current and magnetic filaments is plotted in the z=150 and z=350 subsections. 
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FIG. 8. — Longitudinal heating and acceleration of ambient ions. The top {bottom) panels are for a simulation in which the initial momentum distribution 
function of the injected pairs is monoenergetic (broadband) and the ambient medium is composed of electrons and ions. Also shown is the average transverse 
magnetic field amplitude (in arbitrary units) in the x-y plane as a function of z (solid curves). 
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FIG. 9. — Current filamentation and magnetic field energy density after saturation: t = 59.8u;~ . The amplitude of the average (in the x-y plane) transverse 
magnetic field, B± is plotted (in simulation units) as a function of z together with the average z component of the total current density, J. (top panel). The z 
component of the ion current density is plotted for a subsection of the computational box: y = 43A (middle panel). The a (b) panels are for a simulation in which 
the initial momentum distribution function is monoenergetic (broadband). Color bar gives the amplitude of J z = [/max,— Jmux] and B± = [0,B raax ] in simulation 
Units: Jmnx = 24.6 (22.0), /max = 15.1 (7.8) for ambient ions and Z? max — 3.6 (2.8). 



